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1. Introduction 

While there is no generally agreed definition of gravitational entropy in General 
Relativity (GR), it was conjectured some years ago by Penrose [20], that it should 
be related to the clumping of matter, therefore to the degree of inhomogeneity 
and anisotropy of a space-time and therefore associated with the Weyl or conformal 
curvature. Specifically, Penrose suggested that a measure for the gravitational entropy 
should involve an integral of a quantity derived from the Weyl tensor, and that the 
definition of the scalar product on the Hilbert space of one-particle states for linearised 
GR in flat space might provide guidance [20]. This deflnition, which we review below, 
is written in terms of potentials for the (linearised) Weyl spinor.| 

There have subsequently been several attempts to construct gravitational entropy 
measures using polynomial invariants of the Weyl and Ricci tensors (see e.g. [5, 10, 19]) 
but with no complete success. Prom the point of view of Penrose's suggestion these 
have the wrong 'differential order' being constructed algebraically from the Weyl tensor 
rather than from potentials for it. Another approach has been to note that, via the 
Bianchi identities, the Weyl tensor is related to the density gradient, which is a natural 
physical measure of inhomogeneity. There are proposals for a definition of gravitational 
entropy based on density contrast functions which are covariant, non-local, and work 
globally for dust cosmologies, provided the initial singularity is isotropic [12, 16]. This 
is encouraging, as part of the motivation of [20] was to connect the notion of low 
gravitational entropy to a restriction on the nature of the initial singularity. 

It is unclear how any of these measures could relate to the established notion of black 
hole entropy or whether there should be a relation between the entropy in gravitational 
waves and black hole entropy at all. 

Our purpose in this paper is to return to Penrose's suggestion and take it as far as 
possible, using the scalar product from linear theory in fiat space to motivate a definition 
in curved space. This requires a potential for the Weyl tensor or Weyl spinor and we shall 
use the Lanczos potential ([14], see also [1], [3]). It is a general result of Illgc [13] that 
any spinor field with the symmetries of the Weyl spinor locally has a Lanczos potential 
which is (uniquely) determined by its value at a space-like hypersurface S. Furthermore, 
for a vacuum space-time there exists a potential for the Lanczos potential, a second or 
super-potential for the Weyl spinor, again determined by its value at S [13] (see also 

[I])- 

Apart from lUge's result, there is no general prescription for obtaining a Lanczos 
potential for a given spacetime. A general expression for a Lanczos potential in the case 
of perfect fiuid spacetimes with zero shear and vorticity was given in [18]. More recently, 
this result has been extended by Holgersson [11] to Bianchi I perfect-fiuid spacetimes. 
There are also several examples in the literature of Lanczos potentials for particular 

I For a different approach to implementing Penrose's suggestion via graviton number in FLRW 
cosmologies, see [17]. For more on potentials in linear theory and electromagnetism see e.g. chapter 6 
of [21]. 
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exact solutions, including Godel, Schwarzschild, Taub and Kerr [4,6,7,15,18]. Since 
there is freedom in the choice of the Lanczos potential, these examples are usually made 
subject to symmetry assumptions. 

In this paper, we consider the Lanczos potential and superpotential for linearly 
perturbed Friedman-Lemaitre-Robertson- Walker (FLRW) spacetimes. We obtain wave 
equations for tensors defining invariant parts of the Lanczos potential, and solve them 
by Fourier transforms to give explicit solutions. Then, as an application, we closely 
following Penrose's idea [20], in order to propose a measure of the gravitational entropy 
and apply it to linearly perturbed FLRW cosmologies. The measure is defined as far as 
possible in such a way as to carry over to more general situations. Thus given a choice of 
'time', in the sense of a space-like hypersurface, we construct a Lanczos potential with 
data at that hypersurface, and a tensor which obeys the equations for a second potential 
only at that hypersurface. Then we define a complex structure at the hypersurface on 
the space of potentials. From this, we construct a measure of gravitational entropy Sg 
which is now a functional of the (linearised) Weyl tensor. 

Our definition is of course speculative. It has a reasonably good motivation, but 
doesn't have any obvious monotonicity property. It has at least one good property: 
the main positive result of this part of the paper is that Sg is finite at the initial 
singularity only for those linearised Weyl tensors which are finite at the initial singularity. 
Here the Weyl tensor initially is understood to be finite if the metric (background plus 
perturbation) can be rescaled so as to extend conformally through the singularity, in 
other words if the initial singularity is still isotropic [10]. If the initial Weyl tensor is 
singular, in the sense understood here, then the initial gravitational entropy is infinite. 
Thus finite initial gravitational entropy, as defined here, requires finite initial Weyl 
tensor. 

The plan of the paper is as follows. To end this section we review linearised 
GR, introducing the scalar product on the space of solutions and explain how this 
motivates Penrose's suggestion. In the next section, we review the (1 + 3)-formalism 
for cosmological models and apply it to perturbations of FLRW cosmologies. Then 
in Section 3, we obtain an expression for a Lanczos potential for these perturbations, 
solving the wave equations which arise by Fourier transforms. Finally in Section 4, we 
use the analysis from earlier sections, together with a definition of complex structure 
on the space of potentials, to suggest a definition of gravitational entropy for these 
cosmologies. 

In linearised GR, we start by perturbing the fiat metric rjab according to the equation 

9ab = Vab + ^ab- 



For simplicity, we shall assume that the perturbation ^ab is subject to the following 
gauge conditions: 

= = V"$„5. (1) 
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We obtain the linearised connection (in the sense of the perturbation of the Ricci rotation 
coefficients) as the tensor L^bc — -^[a&]c defined by 

Labc = V[„$6]c. (2) 

Note that 

Lab' = = V'^'^'Labc = V'L^bc (3) 

With the convention 

the hnearised Riemann tensor is 

Rab''' = -'^laL''\]-^^''Lj. (4) 

The field equation is the linearisation of the Einstein vacuum equation and is equivalent 
to the wave equation on ^ab- 



In spinors we write = ^aba'B' (not to be confused with the Ricci spinor) and 

Labc — LabCC'^A'B' + C.C. 

for a symmetric spinor Labcc- Then (2) and (4), taking account of (1), are 

^(f ^ABC'D' — —'^Labcc'i (5) 

'^D Labcc =iPabcd, (6) 

where iPabcd is the Weyl spinor. 

We define the symplectic form on the space of solutions by 

n{^, $) = / {Labc^'^ - labc^'")dS\ (7) 
J S 

where the integral is over a space-like hypersurface S. Then this is independent of surface 
by virtue of the field equation, and is also gauge-invariant. For the scalar product on 
the Hilbert space of classical solutions, we need the complex structure J on the space of 
solutions, which is usually defined by the splitting into positive and negative frequencies: 
if $ = $^ + $_ is that sphtting then 

J$ = i($+-$_), 

and then the inner product is 

0(*, J*). (8) 

This can be seen to be positive-definite by writing it in terms of the Fourier transform. 

This construction motivates Penrose's [20] suggested guide to a definition of 
gravitational entropy: given a solution to a classical field theory, the coherent state built 
on that classical state can be thought of as the most closely corresponding quantum 
state; the expectation value of the number operator in the quantum state is then a 
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measure of the 'number of particles' underlying the classical state, which in turn is a 
measure of the entropy; but this expectation value is just the norm of the classical state 
in the scalar product. If it were possible to find a definition like (8) above in a curved 
space-time, it would therefore be a candidate for a definition of gravitational entropy. 
(Again, this might be just the entropy in gravitational waves. Black hole entropy could 
be another story.) 

This is the idea we pursue here. We replace the linearised connection L^bc by 
the Lanczos potential, which by lUge [13] always exists satisfying (6). The second or 
superpotential derived from Labc always exists for vacuum but we shall be concerned 
with cosmological solutions when its existence is problematic, as is the correct definition 
of J. On the other hand, we do not expect a definition which is independent of surface 
or, equivalently, independent of time. Thus our aim will be to mimic (8) as close as 
possible but basing the construction on a choice of hypersurface. 



2. 1+3 Formalism and tensor perturbations of FLRW cosmologies 

We consider a space-time with a distinguished time-like direction given by the velocity 
vector field u"' of the fluid, and use the formalism of [8,9,22], with gabU°'u'^ = —1- We 
introduce the tensor which, at each point, projects into the space orthogonal to u°' by 

Kb = Qab + UaUb. (9) 

Then 

Khl = hl h':ub = o, K = 3. (10) 

The covariant derivative of can be written, as usual, as 

VaWfe = ^Ohab + (Tab + ^ab " ^Ub (11) 

where 

(Tab = Cr(ab)'-, CTa = 0; (T abu'' = 0; Uab = OJ[ab]'-, ^abU^ = 0. (12) 

Then it" is the acceleration, Uab is the vorticity tensor, aab the shear, and 6 the expansion. 
The stress-energy tensor for perfect fluids is 

Tab ^ pUaUb + phab (13) 

where p is the energy density and p the isotropic pressure of the fluid. 

We shall be principally concerned with the case of vanishing vorticity. Then the 
fluid flow is orthogonal to space-like hypersurfaccs St, which can be labelled by proper- 
time t along the flow, hah is the (Riemannian) metric on these hypersurfaccs and its 
Lcvi-Civita covariant derivative, say Da, is defined by projection: if Vb...c is a tensor 
orthogonal to u'^ on all indices then 

DaVb...c^hihl---hiVdVe...f. 
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A useful operation below will be curl, defined for a symmetric tensor Xab orthogonal to 

u"' by 

(curl Xy'' := r]'^^''D^X%, (14) 

where r]abc — Vabcdu'^ is the volume form of S^, and rjabcd is the space-time volume form. 
It will frequently be convenient to omit the brackets and write just curl X"-^ or curl X^b- 
The Weyl tensor can be decomposed into its electric and magnetic parts, Eab and 
Hab relative to the velocity vector u"' as 

Eab = CacbdU'^u'^, Hab = CacM«"«^ (15) 

where the dual Cacbd is 

Cacbd = -jVafClstbd- (16) 

From their definition, Eab and Hab are symmetric, trace-free and orthogonal to u"". 

The Bianchi identities for the Weyl tensor in the case of a twist-free perfect fluid 
which also has vanishing acceleration can be written as the following system (see e.g. 
[22]): 

+ Za,^''E^^' -(TcdE'^h''\ (17) 

+ 3a^^''H^^'' -acdH'^h''^, (18) 
^ry^'^V/^^rf+^^V, (19) 

DaH^' = - v'''(^ceE^d, (20) 

where the dot is m^Vq. 

Now we use this formalism to consider perturbations of FLRW cosmologies which 
are purely gravitational. The background is conformally-fiat, so that Eab = Hab = 
and the fluid-flow is geodesic, shear-free and twist-free so that 

Ua = UJab = CTab = 0. 

We consider the FLRW metric gab linearly perturbed with Sgab = ^ab- Following [23], 
for purely gravitational perturbations we may consistently impose the gauge conditions 

^abu' = ^\ = V^^ab = 0. (21) 

For the linearised field equation, we characterise the perturbation as purely gravitational 
by requiring that the perturbation in the Ricci tensor in the form i?^ vanish: 

SR'a = 0. (22) 

This implies that Sp — 5p — 0, and with the gauge conditions (21) also 6u"- — 5ua — 0, 
so that 5T^ — for the stress-energy-momentum tensor. 
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For the perturbation in the kinematic quantities it easily follows, for example in 
the formalism of [22], that 

69 = = 5uJab = 5ua 
while, for the shear, we introduce the notation: 

T^ab ■= Saab = ^^ab- (23) 

For the Weyl tensor, which is zero in the background, we find (from the equations in 
e.g. [22]) 

^ab ^ _±ab _ (24) 

3 

H"^ = curl E"^ (25) 
Now the field equation (22) is 

U^ah = \p ^ab, (26) 

(compare e.g. [23]). We note the following identities for trace-free, symmetric tensors 
Xab orthogonal to u"': 

D'^cmlxab =0, (27) 

V^Xafe = D'^Xab. (28) 

(V^Xafe)- = ^"Xab, (29) 

(curl Xah)' = curl Xab - -^0 curl Xa&, (30) 

curl curl Xab = - Xab - ^Xab -0 Xab+{p- ^0'^)Xab- (31) 

Then, from (21), (23), (28) and (29) 

D«E„, = 0, (32) 
and from (23) and (26) we calculate 

aj:ab = letab + {lp-lp + le^)^ab- (33) 

o o 2 3 

It is now easy to check that, neglecting second-order terms, (17)-(20) are satisfied 
with Eab and Hab as in (24) and (25). Specifically (18) follows from (30), (19) from (29) 
and (32), (20) from (27), and finally (17), which is the hardest, from (26), (31), (33) 
and the Raychaudhuri equation: 

+ j + l{p + Sp)^0. (34) 
3. The Lanczos potential 

The Lanczos potential is a tensor Labc — —Lbac: connected to the Weyl tensor by the 
equation: 

C,, = -VHJ - V^aL^,^ - traces, (35) 
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which should be compared with (4) (in general, we follow [11] but our definition of Late 
is twice the usual definition, in order to maintain (2)). There is gauge freedom in Labc 
satisfying (35), which can be reduced by imposing the Lanczos gauge conditions: 

Lab' = = ri^^'^L,,, = VeL„,^ 
the same conditions as in (3). When these are imposed, the '-traces' term in (35) is 
-^S^la^Qb]^ where 

Qac — ^ -^abc) 

which is symmetric and trace-free by virtue of the gauge-conditions on Labc- This term 
vanishes in the Minkowski space version of the theory described in Section 1, but doesn't 
necessarily vanish in curved space. 

The algebraic gauge conditions ensure that Labc can be expressed in terms of a 
symmetric spinor field Labcc ^ 

Labc — LabCC'^A'B' + La'B'C'C^AB, 

and the differential gauge condition then implies 

Labcc = 0. (36) 

Now (35) takes the spinor form 

V^^ Labcc — iI^abcd, (37) 

just as in (6) but where iPabcd is now the full (nonlinear) Weyl spinor. There is no need 
to symmetrise in (37) because of (36). lUge [13] shows that (37) has a unique solution 
given Labcc on a space-like surface, but one cannot in general find the second potential 
as in (5) as this equation has a curvature obstruction from the Ricci tensor: given (6), 
(5) implies 

f^''''%^B)ECD' = 0, (38) 

where ^aba'b' is the Ricci spinor. 

Holgersson [11] gave a useful decomposition of the Lanczos potential into irreducible 
parts in the (l-|-3)-formalism as follows: 

Labc = 2U[aAb]Uc - A[ahb]c - '^■U[aCb]c + Vab'^Sdc + U[arjb]cdP'^ - UcfjabdP'^ ,{^9) 

where Aa and Pa are orthogonal to and Sab and Cab are trace- free, symmetric and 
orthogonal to u^. This gives sixteen components for Labc (three each for A and P; 
five each for S and C) which agrees with the eight complex components for Labcc- 
Holgersson [11] also gave formulae for the electric and magnetic parts of the Weyl tensor 
which are useful below. 

We want to calculate a Lanczos potential for a perturbed FLRW spacetime with 
some given perturbation ^ab in the metric, as considered in Section 2. Since the 
perturbation is characterised by a trace-free, symmetric tensor orthogonal to m", we 
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seek a Lanczos potential as in (39) with vector parts zero: Aa — Pa — 0. Then from 
(39) and (35), or quoting from [11] we find 

Eab^^(cur\Sab-Cab), (40) 

^^{curl Cab + Sab). (41) 

Equating these to (24) and (25) we have equations for Cab and Sab- Now, if a 
superpotential 0a6 also existed for all times with 

Labc = ^[a4>b]c (42) 

then, from (39), we would have 
1 6 

Cab = ^{<i>ab + g^ab), (43) 
Sab = ^CUrl (pab, 

but this is incompatible with the Bianchi identities (17)- (20). This incompatibility is a 
consequence of the obstruction (38). However, (43) suggests another ansatz, namely 
1 9 

Cab = ^{'ipab + g0a&) 

Sab = ^CUrl (f)ab (44) 

in terms of another unknown tensor ipab- (This is simply an ansatz, in that we express 
two unknown tensors, C and S, in terms of two other unknown tensors, (j) and ip] the 
justification of the ansatz is the simplification which results, for example in (50) and 
(51) below). 

Using (30) we find from (41) and (44) 

Hab = |CUrl (0 + V')a6- (45) 

Comparing this equation with (25) we can choose 

^ab^^{<Pab + 1pab), (46) 

so that ipab is known once (pab has been found. Then, using (40) 
1 • 6 6 ■ 

Eab = -^{-'^ab - -^(l>ah " g^afe + CUrl CUrl (pab), (47) 

and combining this with (24), (31) and (46) we get 

4 9 9"^ 8 

^(t>ab + 3^0a& + (3 + -g - P)'t>ah = g^^^b, (48) 

which is a wave equation for (pab- Note that this wave-equation is not (26). In fact, if 
we introduce Xab — (pab — ^ab then Xab satisfies 

4 9 9"^ 1 

nXab + ^9Xab + (3 + - P)^ab = l^iP + P)^ab, (49) 



Lanczos potentials and gravitational entropy for perturbed FLRW space-times 



10 



using the Raychaudhuri equation (34) again. Equation (49) is not satisfied by zero, so 
that (pab cannot be taken to be the perturbed metric For later use, we note that, 
in terms of Xab and 

(f>ab = ^ab + Xab, (50) 
Ab = ^ab - Xab- (51) 

Given initial data (0„6(x, to), 0ab(x, to)) or equivalently (0„6(x, to), V'a6(x, to)), in terms 
of spatial coordinates x at some time to, a solution to (48) exists and is unique. We 
therefore have a complete prescription to determine a unique Labc for linearly perturbed 
FLRW, subject to this data. We can achieve (43), and therefore (42), at a given instant 
by choosing the data to be 

0„6(x,to) = $a6(x,to) (52) 
0a6(x,to) = ^a6(x,to) 

at that instant, or equivalently 

X„5(x,to) =0 = X„5(x,to), (53) 

but this won't then be true at other times. 

We summarize our results so far in the following proposition: 

Proposition: Given a perturbed FLRW spacetime and a choice of time to, a Lanc- 
zos potential Labc, in the Lanczos gauge, may be uniquely specified by (39) with (44), 
(46) and (48), subject to the data (52). We may define a superpotential (pab such that 
(42) holds at to but this will not hold at other times. 

To obtain exphcit solutions for and X, wc first recall some details of the FLRW 
metrics. For simplicity, we shall assume A; = 0, the spatially-flat case, so that the metric 
is 

g = R{tf{dx^ + dy"^ + dz^) - dt^, 

or, introducing conformal time r, 

g = Rirfidx'' + dy'' + dz'' - dr''), 

where R{t) = R{t) and dt/ R(t) = dr. The overdot will consistently stand for d/dt, and 
d/dr will always be written explicitly. 

We choose, as spatial coordinates, x = (x') for i = 1,2,3, and note that dS — R^d^x 
is the volume element on the hypersurfaces of constant t. We shall assume an equation 
of state of the form p = (7 — l)p, and then 9 = 3R/R, while the conservation equation 
implies that p — poR~^^ for constant po- The Friedmann equation reduces to 

if = (54) 

(with the convention that SttG/c^ = 1) and, without loss of generality, the solution is 
R = t" where n — 

37 
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In coordinate components, for a tensor Xab orthogonal to u"- the d'Alembertian is 

1 R R R^ 

= -^^oXij - Xij + -^Xij + 2(-^ + -^)Xir (55) 

where Aq is the fiat Laplacian in x, and for the time-evolution we find, in components: 

Define the Fourier transform l>ij(i,q) of the coordinate components x) of $a6 in 
the usual way as 

then, suppressing indices for clarity, (26) with the aid of (55), becomes 

R ^R R^' ' 
which, as a check, is equation 15.10.39 of [23] (when making the comparison, recall that 

R \, 

^ = -g(p + 3p) 

with our conventions). Following [23], substitute 

^ij = r"(/..^(q)^+(r, q) + h^{c,)E^{T, q)), (57) 

where a. — 2^fi^ and /iy(q) and /lij(q) represent the two polarisation states i.e. they 
are symmetric, trace-free matrices of a standard form, orthogonal to q and suitably 
normalised. Then (57) gives 

+ r— + (lql'-:^)^ = o, 



dr'^ r dr 



where u — ^jz^, which is Bessel's equation of order u in |q|r, so that each of H'^ and 
if ^ is a q-dependent linear combination of Bessel functions J±i/(|q|T). 

Now we need the Fourier transform of (49), in coordinate components, taking 
account of (54), and again suppressing indices for clarity. This is 

5R^ .3R 9i?2 |q|2 i 

Analogously to (57), put 

^ij = r/^(/^J(q)G+(r, q) + ^](q)G^(r, q)), (59) 
with (3 = to find 

^+i^+(|<.r-3i.)G^-^."-^^ (60) 

and similarly for in terms oi . 
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We may solve (60) by variation of parameters: suppose t — to corresponds to r = a, 
then we want G{a) = G{a) = by (53). The homogeneous equation is Bessel's equation 
with V — 1/2, and the solution for is 

^^(^) = Tr^^\ / o-^-'l^ ^x^{H{t - a))H+{a)da, (61) 
— n) |q| Ja 

and similarly for (7^. 

To summarise: (57), (59) and (61) determine X and $ and then (50) and (51) 
determine and ip; finally (44) and (39) determine the Lanczos potential. 

4. Gravitational Entropy for perturbed FLRW 

We want to recapitulate as much as possible of the argument in Section 1, using the 
Lanczos potential in place of the L used there. In place of (7), given two solutions of 
(26), the Einstein equations for gravitational perturbations of FLRW, we shall define 

*(2)) = / {l'^U'-'^"^ - L^L<t>^'^"^) u^dS, (62) 

with (f)^^l related to L^^l^ as in (42) at a particular choice of to- Then from (39) and (44) 
we find 

= 1 1 i^u^'^''' - 0i;v(^)"^) ds. (63) 

Unlike the case of linear theory, this integral is not independent of time, since 
From (63), using (50) and (51) 

l^($(l),$(2)) = i^^($(i)„,$^2') - $(2)a6$S) dS 

- 2 ~ ^(2)a6^ri)) dS. (64) 

Define the current 

then by (26) j| is conserved: 

V„4 = 0. 
Therefore the integral 

/ UajidS = - f ($(l)a6$^2') - (2)ab^'il)) dS (65) 
J St J St 

is constant in time i.e. the first integral in (64) is necessarily constant in time, but the 
others won't be. (It might be argued that one should use just (65) as the definition 
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of the symplectic form, precisely because it is independent of time; however this is a 
definition which would not be available in more general settings, whereas (62) would.) 

Our strategy now is the following: we seek a definition of gravitational entropy 
based on (8); we don't have a symplectic form which is independent of time, nor do we 
expect to obtain a definition of entropy which is independent of time; however we do 
have a unique definition of first and second potential for the Weyl tensor given a choice 
of time; therefore we shall make a natural choice of complex structure, given a choice 
of time, and use (8) at that time to define the entropy. By using (63) precisely at to we 
simplify (64) greatly, since Xab vanishes there. 

Let us for simplicity restrict to the radiation equation of state, so that 7 = 4/3, 
and then n = 1/2, a = 3/2, (3 = 7/2 and p = 1/2, and we can take R = t^/^ = r/2 and 
pR^ = 3/4. The Bessel functions J±i/2 can be written in terms of elementary functions, 
and as solutions for H^"^^ we can take 

^« = T-i/'(^«(q)sin(|q|r) + B«(q)cos(|q|T)). (66) 

Reahty of ^ab implies that 

l>(t,q)=$(t,-q) (67) 

and then from (57) and (66) the same holds for A^'^'> and B^'''\ 

Note that A parametrises a 'growing mode' and B a 'decaying mode' in the standard 
terminology, and note also that a growing mode is a perturbation of the conformal metric 
which is finite at the initial singularity while a decaying mode is not (essentially because 
for the former $ = 0(t) as i — > 0, just like the unperturbed spatial metric, while for 
the latter $ = 0(t^/^), which diverges by comparison with the spatial metric). Thus 
perturbations with nonzero B have infinite initial Weyl tensor while perturbations with 
B zero but A nonzero have finite initial Weyl tensor. This will be important below. 

Introduce 

(q) = (q)B(2) (-q) - 5(1) (q)^(2) (-q) (68) 
then the first integral in (64) is 

I ^/^(i)"^^^ - = / h\K{^)d\ (69) 

which, as expected, is constant in time. 

For the complex structure on the space of solutions of (26), we follow the discussion 
in [2]. There, for Klein-Gordon fields, the symplectic form is taken to be 

Q(((/),7r),(0,7f))= / i7r4>-<f>n)dS, 

J St 

with TT = 0, and the complex structure is defined by 

J(f) = (6)-^/^ ; Jtt = -(6)^/^- 

Here = —DaD"^, where Da is the intrinsic covariant derivative on St- We seek to 
follow this prescription, with (63) as the symplectic form so that ip then takes the role 
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of TT. Also = — i?~^Ao where, as before, Aq is the flat , three-dimensional Laplacian. 
In terms of the Fourier transform, = |qpi?~^ so that J becomes 

J0 = -^V; ; J^^-\fl (70) 
|q| rC 

We use (57) and (66) to translate (70) into an action on A and B, all evaluated at r = a. 

First for the relation of and -ijj to A and B from (57) and (66) we have 

( 4>\^( Mn Mu\(A\ 
) [m2i M22 )[b ) 

where 

Mil = sin a|q| 

M12 = acosa|q| 
2 

M21 = -(sina|q| + a|q| cosa|q|) 

a 

2 

M22 = -(cosalql — alql sin alql) 

a 

Then, conjugating J defined by (70) with M and introducing a :— a\q\, we find the 
action of J on A and B to be 

( J^\^( Jii Ji2\(A\ 

[jB W21 ^22 M 5 / 



where 



11 

Jii = cos 2q; — — r sin 2q;, (71) 



J\2— r(cos q; — asin a)^ — cos^ a, 

J21 = — (sina + a cosa)^ + sin^ a. 



'^22 — — Jii- 
We note the behaviour as a — > 0, when 



Jii= -- + 0{a), (72) 

a 

Jl2= -^ + 0(1), 



a" 

,2\ 



J21 = 4 + 0(q;'). 

In (68), following the analogy of (8), we set A^^) = JA^^) and S^^) = JB^^\ to find 



x(q,a) = j2iA«(q)^«(q) + J22A«(q)5('Kq) 

- JnS(i) (q):4a)(^ - Ji25(i) (q)SW(^. (73) 

where the choice of J has introduced an explicit dependence on a. Substituting this 
into (69) gives our definition of gravitational entropy at proper time to, corresponding 
to conformal time t — a: 

S,{^- to) = J |q|X(q, a)d\, (74) 
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with K{q, a) as in (73). 

We deduce the following properties of our definition: 

• The entropy at time to is determined by the Weyl tensor at to, but not locally (since 
e.g. the operator (0)^^^^ in (70) is non-local). In particular, it is not the integral 
of a scalar invariant of curvature. 

• Prom the specific form (71) of J (though it follows more generally) Sg can be seen 
to be positive definite. 

• However, for small a we note from (72) that if 5 7^ then K — 0(a"^), 
which diverges on the approach to the initial singularity, while ii B — but 
A then K — 0(1), which is finite. In other words, with this definition the 
initial gravitational entropy is finite if the initial Weyl tensor is finite, and infinite 
otherwise. This is an important success for the definition. 

• Necessarily, for a decaying mode the gravitational entropy will decay from its infinite 
initial value and tend to a positive constant in the remote future (for this, from (73), 
we note that J12 in (71) tends to a constant as t ^ 00). But J12 isn't monotonic 
in time, which makes it unlikely that Sg is. 

• From (71), J21 runs from a value of 4 at a = to a limit of 1 at large a and 
is also not monotonic (though it is positive). Thus for a purely growing mode, 
the gravitational entropy runs from a finite initial value to a finite but lower final 
value, and again there is no reason to expect it to be monotonic. This is perhaps 
less successful for the definition, but this is after all a result from linear theory, 
though in a non-flat background - without nonlinearity, there is no gravitational 
clumping. 

The deflnition has been framed in such a way that it should extend to more general 
situations, and it remains to be seen whether the property of being initially flnite only 
for an isotropic singularity persists. 
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